Hanbury Brown-Twiss Interferometry for Fractional and Integer Mott Phases 
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Hanbury-Brown-Twiss interferometry (HBTI) is used to study integer and fractionally filled Mott Insulator 
(MI) phases in period-2 optical superlattices. In contrast to the quasimomentum distribution, this second order 
interferometry pattern exhibits high contrast fringes in the insulating phases. Our detailed study of HBTI sug- 
gests that this interference pattern signals the various superfluid-insulator transitions and therefore can be used 
as a practical method to determine the phase diagram of the system. We find that in the presence of a confining 
potential the insulating phases become robust as they exist for a finite range of atom numbers. Furthermore, we 
show that in the trapped case the HBTI interferogram signals the formation of the MI domains and probes the 
shell structure of the system. 

PACS numbers: 71.30.+h, 03.75.Lm, 42.50.Lc 
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I. INTRODUCTION 

Cold atoms in optical lattices are becoming ideal many- 
body systems to attain laboratory demonstrations of model 
quantum Hamiltonians due in part to the dynamical experi- 
mental control of the various parameters at a level unavail- 
able in more traditional condensed matter systems. Recent 
experiments with bosonic atoms have been able to simulate 
effective one dimensional systems nl Ij. B "SI, to enter 
the strongly correlated Tonks Girardeau (TG) regime |2, 3] 
predicted many years ago |a] and to realize the superfluid to 
Mott insulator(MI) fT'] transition by tuning the lattice param- 
eters 1 8]. In fact, one of the most active frontier in cold atom 
systems is to explore the possibility of creating new quan- 
tum phases|9]. The capability to physically create them de- 
mands sophisticated diagnostic tools for their characterization 
and actual observation in the laboratory. In this paper, we 
use second order interference techniques, namely the Hanbury 
Brown and Twiss interferometrv llOll to describe the rich phase 
diagram of interacting bosonic atoms in the presence of two 
competing lattices. 

Multiple well super-lattices have been experimentally real- 
ized by superimposing two independent optical lattices with 
different periodicities] 11, 12]. Recent theoretical studies of 
the many-body Bose-Hubbard (BH) Hamiltonian in the pres- 
ence of an additional superlattice 1 13 nl4il5lfl6lfl7lflal have 
centered mostly on the detailed phase diagrams of the system. 
The focus of this paper is to characterize the different phases 
of the system by analyzing the four and two point correlations 
that can be extracted from the time of flight images. Our study 
is confined to period-2 superlattices, the simplest example of a 
system exhibiting almost all the key aspects of a more general 
period-TTT. lattices. 

The HBTI technique measures second order correlations 
from the shot-noise density fluctuations in the absorption im- 
ages of exp andin g atomic gas clouds. This technique has been 
shown rl^ l20ll2lll to be particularly suited for probing many- 
body states of cold atomic systems, as it provides complemen- 
tary information to the first-order correlations inferred from 
the average density distribution in the absorption images. For 



example, cold bosonic atoms in the MI regime do not exhibit 
first order interference pattern but they do have sharp second- 
order Bragg peaks which reflect the spatial periodicity of the 
lattice. The interference peaks are the manifestation of the en- 
hanced probability for simultaneous detection of two bosons 
(bunching) due to the Bose-Einstein statistics. 

The phase diagram of interacting bosons in the presence 
of ID superlattices is landscaped by various quantum phases. 
Transitions to the different phases can be driven either by 
changing the ratio between the three sets of energy scales 
in the system: the interaction energy U, the tunneling rate 
J and the the lattice modulation A or by varying the filling 
factor V. The lattice modulation induces additional fractional 
MI phases that occur at filling factors commensurate with the 
periodicity of the combined superlattice. These phases can 
be understood in the hard core boson (HCB) limit where the 
strongly interacting bosons can be mapped to noninteracting 
fermions. The superlattice fragments the single particle spec- 
trum, inducing band gaps. The filling of a sub-band at a crit- 
ical filling factor in the fermionic system results in a band in- 
sulator state that corresponds to a fractionally filled MI state 
of the bosonic atoms. The fractional MI phases survive in the 
soft core boson limit beyond a critical value of the interaction 
energy U. 

Besides the fractional Mott phases, the interplay between 
the interaction energy and the superlattice potential can also 
lead to a new type of integer-filled Mott phases which cannot 
be understood in any special limit such as the HCB limit or 
the pure periodic case as they appear at finite U values and for 
finite strengths of the superlattice potential. These insulating 
phases are characterized by a modulated density profile and 
hence will be referred as staggered Mott phases. 

Here we demonstrate that HBTI provides definite means to 
monitor the various superfluid to MI transitions and we use it 
as a method to obtain the phase diagram. Firstly, for a fixed 
U, J and A we vary f, and find that the onset to the transi- 
tion is accompanied by a change in the sign of the superlattice 
induced Bragg peaks as the filling factor is increased beyond 
the critical ones. Secondly, for a fixed filling factor ly and fixed 
A/ J, we change the ratio U/J across the transitions and show 



that they are signaled by a sharp maximum in the intensities 
of the Bragg interference peaks. We show that the phase dia- 
gram obtained from HBTI is in qualitative agreement with the 
recently reported phase diagram calculated from Monte-Carlo 
simulations 1 14, 15. 16. 17. 181 and with a mean field phase 
diagram that we derive analytically. 

Another important aspect of our study is to investigate the 
effects of a harmonic trap on the various phases. In contrast 
to the translationally invariant system where the transition to 
insulating phases occurs only at critical filling factors, the 
parabolic confinement allows for the formation of fractional 
and integer filled MI domains in a finite window of fillings. 
We show that harmonically confined superlattices in the HCB 
limit develop a shell structure |22] analogous to the one ob- 
served in the absence of any lattice modulation at moderate 
values of the interaction. The HBTI interferogram provides 
clear signature of the formation of the different Mott domains. 

The paper is organized as follows. In Sec.|II|we define noise 
correlations and the model Hamiltonian used in our study. In 
Sec. Hn] we calculate two and four-point correlations in the 
HCB case. In Sec lIVI we consider the BH system for finite 
values of interaction. We first analytically calculate the phase 
diagram by using mean field theory ( whose details are de- 
scribed in Appendix A) and compare it with the phase diagram 
obtained from noise correlations. These higher order correla- 
tions are calculated by numerical diagonalization of the BH 
Hamiltonian. In Sec|V] we discuss the effects of the harmonic 
trap in the HCB limit. In Sec lVII we state our conclusions. 



II. NOISE CORRELATIONS 

In a typical experiment, atoms are released by turning off 
the external potentials. The atomic cloud expands, and is pho- 
tographed after it enters the ballistic regime. Assuming that 
the atoms are noninteracting from the time of release, proper- 
ties of the initial state can be inferred from the spatial images 
I13l U^ 1^]: the column density distribution image reflects 
the initial quasimomentum distribution, n{Q), and the density 
fluctuations, namely the noise correlations, reflect the quasi- 
momentum fluctuations, A{Q, Q'), 



i{Q) = l^e'«-(^-'=)a] 



a-k, 



(1) 
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A{Q,Q') = {n{Q)n{Q')) - {h{Q)){h{Q')). (2) 

In Eq. (|2} we have assumed that both Q, Q' lie inside the first 
Brillouin zone. Here L is the number of lattice sites and a 
the lattice constant. In this paper, for simplicity, we focus on 
the quantity A(Q, 0) = A((5). In our discussion below, N 
denotes the number of particles in the system and v — N/L 
is the filling factor In this paper, we will confine ourselves to 
the case with v <1. 

The BH Hamiltonian describes bosons in optical lattices 
when the lattice is loaded in such a way that only the lowest 
vibrational level of each lattice site is occupied and tunneling 
occurs only between nearest-neighbor sites |22]. The ID BH 



Hamiltonian in the presence of a period 2-superlattice is given 
by 



+ X 2A cos{tt j)nj + ^ flffij. (3) 



Here a^ is the bosonic annihilation operator of a particle at 
site j, hj — aJ-aj, and the sum (i,j) is over nearest neighbors. 
The hopping parameter J, and the on-site interaction energy 
U are functions of the lattice depth. The cosine term describes 
the onsite potential generated by the additional lattice with 
twice the periodicity of the main lattice. Here the main lat- 
tice creates the tight-binding system and the second lattice is 
assumed to be a weak perturbation. In this case atoms in the 
super-lattice have a hopping parameter which is independent 
of the lattice position. The parameter A is almost proportional 
to the depth (in recoil units of the main lattice) of the addi- 
tional lattice (see Ref. Illll for details). The last term takes 
into account the parabolic potential with J7 proportional to the 
parabolic trapping frequency. 



ni. HARD CORE LIMIT 

In the strongly correlated regime, when U ^ oo, the BH 
Hamiltonian can be replaced by the HCB Hamiltonian |25J, 
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X 2A cos(7rj)nj + ^ nf 



(4) 



Here bj is the annihilation operator at the lattice site j which 
satisfies [hi, bA = 6ij, and the on-site conditions M = bj^ = 
0, which suppress multiple occupancy of lattice sites. 

HCB operators can be linked to spin-1/2 operators by 
means of the Holstein-Primakoff transformation |24|, which 
maps bosonic operators into spin operators. The Holstein- 
Primakoff transformation maps the HCB Hamiltonian into the 
XY spin- 1/2 Hamiltonian. The latter can in turn be mapped 
onto a spinless fermion Hamiltonian by means of the Jordan- 
Wigner transformation 11251 12611 . 

The Bose-Fermi correspondence can be used to calculate 
various many-body observables of the strongly interacting 
bosonic system in terms of the ideal fermionic two-point func- 
tions which can be written as gim = "^k^o ^i V'm . Here 
V^l are the amplitudes of the single atom eigenfunctions at 
site j and _E(j.) are the single atom eigenenergies: 
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J{rj+i+^--i) + 2Acos(7rj)V]'^ + nj'i,';> = E(^k)i^^ 
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The single -particle spectrum of period-2 superlattice can be 
obtained by decimating every other site of the tight binding 



Eq. (|5|i. The renormalized system at even or odd sites is de- 
scribed by 
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J\i^T+2 + V'r^) + (2 J^ + 4A^)^f ^ = Et,)i^^ 
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(6) 



Assuming periodic boundary conditions, the eigenenergies 
are then given by 



E(k) — ±2^ 



/j^cos(M 



A2, 



(7) 



where k — 0, 1, ... , L/2. The effect of the additional lattice 
potential is to split the band into two different sub-bands, each 
with band-width 2 yJ^ + A^— 2A, separated by an energy gap 
of4A. 

Local observables such as the density distribution and ener- 
gies are identical for the HCB and non-interacting fermionic 
systems. For example, the HCB ground state energy corre- 
sponds to the sum of the first N single-particle eigenstates. 
On the other hand, fermions and HCBs posses differ non-local 
correlation functions. A general formulation to calculate HCB 
two-point correlations has been developed by Lieb and Mattis 
f2^. In our earlier studies 1 23], we have generalized Lieb and 
Mattis formalism and have obtained explicit formulas, involv- 
ing multiple Toplitz-like determinants, to compute the four- 
point correlation functions required to calculate noise correla- 
tions in HCBs. The evaluation of these determinants, whose 
order scales with the size of the system, is in general compli- 
cated and therefore an analytical treatment is difficult. Below 
we will describe the results obtained by numerical computa- 
tion of the explicit formulas discussed in our earlier study. For 
information about the exact formulas and various other rele- 
vant details we refer the readers to our earlier paper I23ll . 
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FIG. 1 : Quasimomentum distribution as a function of the filling fac- 
tor. The calculations were done for a system with A/ J = 1, L = 80 
and boxed-like boundary conditions 

Figs. [2 and |3 show the quasimomentum distribution and 
the noise correlations for all filling factors, < i^ < L The 
quasimomentum distribution exhibits interference peaks at the 
reciprocal lattice vectors of the combined superlattice: a large 
peak at Q = and somewhat weaker peaks at Qa — ±7r, 
induced by the period two modulation. The peaks are a mani- 
festation of the quasi-long-range coherence of the system. At 




FIG. 2: Noise interference pattern as a function of the filling factor 
for a system with A/J = 1 , L = 80 and boxed-like boundary 
conditions. In the inset we truncated the central peak to make the 
negative background and the peak to dip transition more visible 



J/ = 1/2 and 1 the first order interference pattern flattens out, 
signaling the insulating character of the system at these criti- 
cal fillings. 

The peaks at Qa ~ and ±7r also exist in the noise inter- 
ference pattern (Fig. 2), where they are narrower and are ac- 
companied by adjacent satellite dips, immersed in a negative 
background. These satellite dips are the signatures of the long 
range coherence in the second order pattern as they disappear 
in the insulating phases (see detailed discussion in Ref.[23]). 
The satellite dips are clearly seen in the inset. In contrast to 
the quasimomentum distribution, the peaks at reciprocal lat- 
tice vectors continue to exist at the critical filling factors in 
the HBTI pattern. However, their intensity is strongly reduced 
compared to the intensity at other filling factors. 

To highlight the signatures of the insulating phases in the 
noise correlations, we plot in Fig.|3]the visibility of the cen- 
tral and superlattice induced peaks. We define the visibil- 
ity as the intensity of the second order Bragg fringes (noise- 
correlations) normalized with respect to the quasimomentum 
distribution: 



v(g) 



A(Q,0) 
7i(Q)n(0) 



(8) 



The visibility V((3) is a relevant experimental quantity as the 
normalization procedure filters some of the technical noise in- 
troduced during the measurements. 

In the visibility figure one observes the development of very 
sharp peaks at the critical fillings. In addition, an interesting 
definitive signal of the fractional insulating phase is a change 
in the sign of the intensity of the superlattice induced peaks: 
as the filling factor is increased beyond half filling, the peak 
at Qa = ±7r becomes a dip. This effect, which we will refer 
to as peak to dip transition, has its origin in the occupation of 
the second band with few atoms, once the first band is fully 
occupied. In fact, it can be interpreted as a manifestation of 
the fermionization of HCBs. The peak to dip transition is a 




negative fringes survive even for more than one atom in the 
second band. 

Eg. din also shows that the absence or the presence of an 
interference peak at Qa = tt can be used to distinguish the 
fractional and Mott insulating phases: the peak at Qa — tt 
disappears when the system is a unit filled insulator but con- 
tinues to exist at aQ = tt when the system is a half-filled 
insulator On the other hand, for both insulating phases the 
intensity of the central peak is i^(i/ + 1). 



IV. SOFT CORE BOSONS 



FIG. 3: Visibility (see text for definition) of the Q = (dotted line) 
and Qa — tt (solid lines) peaks for a HCB system with X/J = 1 and 
L = 80. 



generic feature of strongly correlated bosons confined in a su- 
perlattice. It is seen for all values of A and, as we discuss 
later, this signature also accompanies the onset to fractional 
Mott state for finite U, i.e. in the soft core boson case. 

A further insight for this peak to dip transition can be gained 
by looking at the noise correlations as A — > oo . In this 
limit, the correlations can be calculated analytically for v — \, 
V = 1/2 and also for the case of one extra atom beyond half 
filling (N = L/2 + 1). At unit filling the ground state in the 
HCB limit corresponds to a unit filled Fock state. In the half- 
filled case, to a good approximation, the ground state can be 
assumed to be a 1/2-filled Fock state l^g) == |1010 . . . 10). 
This is due to the reduced number fluctuations exhibited in 
this limit, as the tunneling between every two sites scales 
as J^/(4A). At filling factor N = L/2 + 1 the ground 
state can be approximated by a 1/2-filled Fock state with an 
extra delocalized particle at the different unoccupied wells. 



with 



these ansatzs, it can be shown that 
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^{Q)n=L/2+1 = 


-j^ + ^^^ha-Sgn), 


(10) 


A(Q)w=L = 


-- + 2dQa. 


(11) 



As shown by these equations, just beyond half-filling, the Q — 
and Qa = tt fringes have opposite sign. The sign difference 
can be understood by calculating A{Q) for a single particle 
(N = 1) or a single hole {N = L — 1). In this simple case, 
noise correlations can be calculated explicitly for all values of 
A: 



A(Q,0) 



Af=l=L-l 



1/ A^ 
4VA2 + 1 



(Sqq - 6qt,) 



(12) 



The explicit formulas for A{Q, 0)n=l/2+i ^nd A(Q, 0)Ar=i 
show on one hand the similarity between these two cases and 
on the other that the filled sub-band does have an effect on 
the extra particle. Actually as can be seen in Fig.|2]and|3j the 



In this section we relax the hard-core constraint and ex- 
plore the interplay between the finite interaction effects and 
the competing periodicity induced by the superlattice poten- 
tial. The key questions that we address are: (1) how a fi- 
nite value of U affects the various phases observed in the 
HCB limit, (2) how generic are the characteristic signals of 
the phase transitions observed in the HCB noise correlations 
as U becomes finite and (3) what novel characteristics of the 
interference pattern emerge as we explore the ( J, f/, A, i^) pa- 
rameter space. 

Before we describe the details of HBTI, we will first use 
mean field theory to gain some insight about the phase dia- 
gram as {U/J, X/J, ly) vary. Although previous studies based 
on numerical Monte carlo simulations have calculated the var- 
ious phases in the soft core regime lMrrxri6lfl7Lil8il . here 
we obtain the phase diagram analytically by using second or- 
der perturbation theory. It is well known that the mean-field 
calculations provide a good characterization of the phase di- 
agram and gives qualitative prediction of the different phase 
transitions even though the exact transition thresholds may not 
be quite correct. We would also like to point out that earlier 
mean field studies lll4l flSl flfl uM of the phase diagram were 
done in a different case, namely when the hopping parameter 
is site dependent. 



A. Mean-Field Phase Diagram 

To calculate the phase diagram we use the well know 
Gutzwiller approximation 1 27 ] that decouples the kinetic en- 
ergy term of the Hamiltonian by introducing a superfluid or- 
der parameter. The details of the mean field calculations are 
discussed in Appendix A. Mean field theory predicts three in- 
sulating phases in addition to the superfluid phase. The insu- 
lating phases can be classified as: (1) a MI phase with i/ = 1 
where all sites have unit occupancy, (2) a staggered MI phase 
with v = 1 characterized by every alternate site doubly oc- 
cupied and (3) a fractional MI state with v = 1/2 where al- 
ternate sites are singly occupied. At other filling factors, the 
system is always a superfluid. 

Fig-Ssummarizes the mean-field calculations. The top left 
panel shows the critical values of U = U/J and X = X/J 
for the onset to the three insulating phases. They correspond 
to the tips of the insulating lobes ( multicritical points) that 
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FIG. 4: (color online) Mean field phase diagram for the period two superlattice. The top left panel shows the critical curves, Ud J — Ac/ J for 
the onset to various transitions to insulating phases. The dash dotted blue line is for the transition to the band half-filled MI The the long dashed 
and short dashed asymptotes corresponds to 2.3/(A/J — 1) and 6.1/(A/J) respectively. The red dotted line corresponds to the transition to 
the unit filled MI and the black solid line to the staggered unit filled insulator. The grid lines are at A = J and 3.75. The top right, bottom left 
and bottom right panels show the phase diagram as a function of the chemical potential for systems with A = 0.5 J, 4 J and 8 J 



are shown for selected values of A in the various panels of the 
figure. 

One of the interesting results of the mean-field theory is 
the prediction of a reentrance to the superfluid phase from the 
staggered MI phase. In contrast, to the critical curves describ- 
ing the threshold to the fractional MI and integer MI phases, 
which vary monotonically in the (C/, A) space, the curve for 
the staggered MI at i^ = 1 bends (see left panel in Fig. |5) and 
becomes a double-valued function. As a consequence, for a 
fixed A > Ac, with the threshold value being approximately 
Ac = 3.75, as we increase the interaction V , the system goes 
from superfluid to staggered MI, reenters the superfluid phase 
and ends in the unit filled MI phase. Note that this reentrance 
to the superfluid phase is absent for A < Ac, where one only 
sees a direct transition from superfluid to unit MI. This behav- 
ior is explicitly shown in the other panels of the figure where 
the chemical potential is plotted as a function of the interac- 
tion energy, for a few selected values of A. In the top right 
panel A < Ac and so only the unit Mott and fractional Mott 
phases are present. For A > Ac, the staggered phase appears 
as an isolated island in the /i — [/ plane whose size grows with 
A. 

The numerical Monte-Carlo calculations reported in Ref. 
1 18] are in agreement with the intermediate superfluid phase 
found at mean-field level. However, according with their cal- 



culations the superfluid region is only present for [7 < 12. For 
[7 > 12 the system is either a unit filled MI or a staggered MI. 
Our mean field calculations, on the other hand, predict that 
the intermediate superfluid phase exists for all values of I] . 
As we will discuss later, this result is in qualitative agreement 
with our numerical calculations performed by exact diagonal- 
ization of the BH Hamiltonian for finite size systems. 

The staggered phase describes an interesting manifestation 
of the interplay between the onsite interaction energy, C/, and 
the energy modulation introduced by the superlattice poten- 
tial, A. This phase is absent in the HCB limit and only exists 
for values of A = A/ J beyond a threshold value and for a 
finite window of [7 = U j J values. The analytic calculations 
predict that the threshold value is Ac = 3.75, which is an over- 
estimation of the corresponding value, Ac ~ 1.75, found with 
the Monte Carlo simulations 1 1^. 

For the half filled case, the mean-field critical curve has 1/A 
as an asymptote (see top left panel in Fig. |3i. This scaling is 
consistent with Eq. (|5) (for large A the effective tunneling rate 
between consecutive odd sites goes like J^/A). Therefore, for 
large A, a fractional MI state can exist for rather small values 
of U . On the other hand Fig. |4]shows that for moderate values 
of A the transition to an insulating phase requires very large 
interactions. We would like to point out that the existence of 
a fractional Mott phase only for values of A > 1 seems to be 



an artifact of the mean-field theory. In fact, in the HCB limit 
([/ -^ oo) the system is a band insulator for any infinitesimal 
value of A. 



B. Noise Correlations and Phase Diagram 

To estabUsh a correspondence between the various phase 
transitions and their signatures in noise-spectroscopy, we re- 
sort to exact diagonalization procedures. We first numerically 
diagonalize the BH system to obtain the ground state of the 
system. This is then used to obtain two and four point corre- 
lations and their Fourier transform. Here we will describe our 
results for 8 wells (L ~ 8). 




FIG. 5: For fixed A/J = 1 and U/J = 20, we plot the visibiUty 
of the central (dotted-dashed line) and superlattice induced (dashed 
lines) peaks as the filling factor varies. The solid and dotted lines are 
the correspondent HCB curves 

Fig.|5ldescribes the visibility [Eq. (|8}] of the second order 
fringes for various filling factors. The important point to be 
noticed here is that the peak to dip transition, a key signature 
of the fractional Mott transition in the HCB limit, is preserved 
for finite U. The figure also shows the corresponding result for 
the HCB case. A comparison between the L = 80 ( Fig.|3} 
and L = 8 results for HCB illustrates the finite size effect on 
the peak to dip transition. As expected, the finite size of the 
system broadens the transition and hence reduces the visibility 
of the interference peaks. Fig. |6] further demonstrates the be- 
havior of the peak to dip effect as the interaction U is varied. 
In Fig|6lwe plot the normalized HBTI pattern as a function 
of U for a system just above half filling. For large but finite 
U > 10, the superlattice induced peak in A{Qa = tt) changes 
sign and becomes a dip. Our numerical analysis confirms that 
the peak to dip transition is a generic signature of the strongly 
interacting regime where bosons exhibit fermion-type charac- 
teristics. 

We next show that the second order interference pattern 
not only complements the characteristics of the various phase 
transitions as observed in the first order Bragg spectroscopy, 
but it also provides new definitive signatures of various 
phases. For clarity, we will describe our results for z^ = 1 and 



Visibility 




FIG. 6: Visibility plot (see text for definition) as a function of U/J 
for a system just beyond half filling (L = 8, A'' = 5) and for A = 2. 



ly ~ 1/2 separately. 
(B.I) Unit-Filled Case 

Fig. Q shows the variation in the central as well as in the 
superlattice induced peak for both the first and the second or- 
der interference patterns as the on-site repulsive interaction is 
varied. An important aspect of the figure is the appearance 
of various local maxima in A(0). The single maximum ob- 
served for small A splits as A increases. A comparative study 
between A^ = 6 and 8 shows that the height at the various 
maxima increase with A^, suggesting that it may be divergent 
in the thermodynamic limit. This observation along with the 
fact that the peak for small A occurs at a value of U very close 
to the well known MI transition, suggests that the maxima 
may be associated with the onset of the various superfluid to 
MI phase transitions. 

By monitoring the locations of the maxima, we obtain the 
phase diagram in the {U, A) plane. Our results are shown in 
Fig- 13 The rather remarkable qualitative agreement between 
this phase diagram and that obtained by mean-field and also 
by earlier Monte carlo studies, supports the validity of our 
conjecture regarding the relationship between a peak in A(0) 
and the onset of a phase transition, and suggests that noise 
correlations can be used as a practical tool to obtain phase di- 
agrams of many body quantum systems. It should be noted 
(Fig. Q that in contrast to A(0), the corresponding zero 
quasimomentum component does not provide any sharp sig- 
natures of the different critical points as U is varied. 

As shown in the figure for A = 3, for moderate on-site 
interaction, the system is in the superfluid phase. As U in- 
creases, correlations begin to build up and when U becomes 
comparable to J, for A > 1.5, the system enters the staggered 
Mott insulating phase with two atoms in the low energy wells 
and none in the high energy ones. The onset to this transition 
is signaled by a peak in A(0) (as seen in the figure) and by 
the development of a positive intensity peak at Qa = tt. It 
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FIG. 7: The figure shows the noise correlations (solid line) and 
quasimomentum fridges (dashed line) for various A values. The left 
and right panel respectively describe the intensities of the central 
peak {Q — 0) and the superlattice induced peak {Qa = n). 
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FIG. 8: Phase diagram for the period two superlattice calculated by 
reading the maxima of the central peak in the A(0) vs. U plot for 
different A . The dotted, dashed and solid lines show the critical 
values in the (7 vs A plane for the transition to the half-filled MI, the 
unit filled MI and staggered filled insulators respectively. 



should be noted that the critical value of A — 1.5 predicted by 
this method for the appearance of the staggered phase is very 
close to the value of 1 .75 found by the numerical Monte-Carlo 
calculations llSll . 

The system continues to be in the staggered insulating 
phase as U increases until the point when the competition be- 
tween U and the energy offset induced by the superlattice, 
4A, drives the system back to a superfluid phase. A simple 



understanding of this reentrance can be obtained by realiz- 
ing that when U is of the order of 4 A the state |2020 . . . ) is 
degenerated with the states |20112011 . . .), 1 10211021 . . .), 
121102110...), and |11201120 . . .) which have superfluid 
character as their average densities are 3/2 and 1/2 in the low 
and high energy wells respectively. The reentrance to the su- 
perfluid phase is signaled by the second peak in A(0). Fur- 
thermore, it is accompanied by the disappearance of superlat- 
tice induced peak at Qa = tt. For values of A close to 1 .5 the 
staggered phase exist only for a very narrow range of U values 
and the first and second order peak can not be resolved. 

As U increases beyond 4A, it is energetically costly to have 
two atoms in the same well and hence the system enters the 
MI phase where sites are singly occupied. This transition is 
also signaled by the third peak in A(0). As A is increased 
from 1.5 up to A « 5, the separation between the last two 
maxima (which determines the range of U values where the 
intermediate superfluid phase exists) decreases. Beyond this 
value, as A is increased further the position of the two peaks is 
shifted to larger values of U. However, the relative separation 
between the peaks was found to remain constant. This finding 
is in disagreement with the results of Monte Carlo simulations 
which find points in C7 — A space where staggered and Mott 
phases coexist. One may argue that the absence of coexis- 
tence between the staggered and Mott phase in our analysis 
may be due to the finite size of our system. Nevertheless, the 
coexistence of such phases seems to imply the existence of 
a first order transition. This makes such a coexistence point 
somewhat subtle and needs further investigation. Besides this 
difference, the phase diagram calculated by monitoring the 
maxima in A(0) is in very good qualitative agreement with 
the Monte carlo as well as the mean field calculations. The 
quantitative differences are due to finite size effects. 

We would like to point out that A{Q = tt/o) also con- 
tains valuable information about the different phases. For 
example the existence of a staggered phase is indicated by 
the development of a positive interference peak at Qa = n. 
This peak is a second order effect and it is not present in 
the quasimomentum distribution. The quasimomentum 
distribution does not give definite signatures of the transition 
points but it gives information about the phase coherence of 
the various phases. As shown in Fig.0 inside the superfluid 
phases the n(0) vs U develops a maximum which disappears 
in the insulating phases where the curve tends to become flat. 

(B.II) Half filled Case 

We next discuss the i' = 1/2 case, where one sees a transi- 
tion from superfluid to fractional Mott phase. The formation 
of this phase is also signaled by the development of a sharp 
peak in the intensity of the central noise correlation peak as U 
is varied across the transition. This is shown in Fig. |9] In our 
finite size system, the peak was found to exist only for values 
of A > 0.205. However, we expect this value to decrease as 
the size of the system increases. Fig. |8]shows the position of 
the peak as a function of A. For large A the critical U value 
decreases as 1/A in consistence with the mean field results. 

The formation of the half filled band insulator is also 
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FIG. 9: The figure shows the noise correlations (solid line) and 
quasimomentum distribution (dashed line) for i/ — 1/2 and various 
A values. The left and right panel respectively describe the inten- 
sities of the central peak Q — and the superlattice induced peak 

{Qa = Tv). 



indicated by the superlattice induced peak at A(7r/a). It 
should be noted that A(7r/a) has negative intensity for the 
non-interacting system. Precisely at [/ = A(7r/a) ~ 
—N/A\^/{\^ + J^). Its amplitude increases as U is increased 
and becomes positive for the values of A that allow the forma- 
tion of a band-insulator. For U > Uc~ , with (Jc~ the 
Mott band insulator critical point, the peak intensity almost 
reaches its hard-core value. 



V. TRAPPED SYSTEM 

In this section we study the experimentally relevant case of 
a superlattice in the presence of an additional parabolic po- 
tential. The two key results of this study are following: First, 
to find that in contrast to the homogeneous system where the 
MI phases only exist at specific critical fillings, the trapping 
potential stabilizes the MI phases and allows for their exis- 
tence over a large range of filling factors. Second, to show 
that noise-spectroscopy provides a detailed information of the 
formation of the insulating domains and the shell structure that 
arises in presence of a parabolic confinement. 

Understanding the effects of the trap on period-2 superlat- 
tice follows along the lines of the earlier studies of the trapped 
system in the absence of any additional lattice modulation. 
Therefore, we first briefly review the key aspects of the A = 
case which has been studied in great detail |28, 29, 30, 31]. 
The combined lattice plus harmonic confinement system pos- 



sesses two distinct classes of eigenstates: low energy states 
that extend symmetrically around the trap center and high en- 
ergy states that are localized on the sides of the potential. The 
origin of these two distinct classes is related to the two energy 
scales in the system, namely the tunneling (J) and the trapping 
energy (il). Modes with excitation energy below 4 J, which 
correspond to the band width of the translationally invariant 
system, are extended and can be thought of as har monic os- 
cillator like modes with effective frequency uj* ~ VUTi and 
effective mass m* — fi,/(2Ja^). Modes with excitation ener- 
gies above 4 J are close to position eigenstates since for these 
states the kinetic energy required for an atom to hop from one 
site to the next becomes insufficient to overcome the potential 
energy cost. The high energy eigenstates are almost two-fold 
degenerate with energy spacing mostly determined by fl. The 
localization of these modes can be understood by means of a 
semiclassical analysis 12911 . Within WKB scheme, the local- 
ization of the higher energy modes can be linked to the appear- 
ance of new turning points related to umklapp processes. In 
contrast to the turning points of the classical harmonic oscil- 
lator that appear at zero quasimomentum, the Bragg turning 
points emerge when the quasimomentum reaches the end of 
the Brillouin zone and can therefore be associated with Bragg 
scattering induced by the lattice. 

The period-2 superlattice splits a band into two subbands 
and hence the nature of modes can be now classified by the 
energy scales associated with these two bands and the band 
gap which is equal to 4A. 

In the top panel of Fig[roiwe plot the various wave func- 
tions. Each eigenstate has been offset along the y axis by its 
energy in units of VjT2. In the lower panel we also show the 
energy spectrum. As shown in the figure, we now see two sets 
of extended and also two sets of localized states. Again, the 
localized states are related to umklapp processes and to em- 
phasize this idea, the classical and the Bragg turning points 
of both bands are explicitly displayed. In Appendix B, we 
provide additional details of the trapped superlattice model. 

Figs. 111! and II 21 describe the density and number fluctua- 
tions that reflect the band structure of the trapped superlattice 
system discussed above. To simplify the description, we in- 
troduce three different quantities: N^^ , Ni and A^^^ which 
respectively denote the number of low energy extended states 
in the first band, the number of states below the first mode of 
the second band and the total number of modes below the first 
localized mode of the second band. In Fig. ^| these three 
number are explicitly indicated with grid lines. 

For filling factors below N^^ atoms tend to spread over 
the central sites populating mostly the low energy wells. The 
extended character of the modes induce large number fluctu- 
ations. For filling factors between N^^ < N < Ni, the 
localized modes in the first band become occupied. Because 
these modes are localized at the edges of the cloud, as the 
number of atoms is increased the occupation of the central 
site remains constant and instead sites farther away from the 
trap center become populated. Thus the presence of single- 
particle localized modes in the Fermi sea leads to the forma- 
tion of fractionally filled insulating domains at the trap center 
in the many-body system. Alternatively, as soon as localized 
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FIG. 10: (color online) The top panel shows the single-particle eigen- 
states of a system with X/J — 1.1 and il/J = 0.005. The height 
of each mode is proportional to its energy in units of VflJ. The 
parabola correspond to the classical and Bragg turning points that 
appear in the first and second band respectively (see text). In the 
lower panel we show the spectrum of the modes in units of \/0,J. 
The dots, triangles, boxes and stars indicate the low energy modes 
in the first band, the high energy localized modes in the first band, 
the low energy modes in the second band and high energy localized 
modes in the second respectively. Note the degeneracy that appears 
between the localized high energy modes of the two bands. The grid 
lines are at TVf ^ = 9, iVi = 32 and TVa^^ = 42. 
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the shell structure observed in the A = case at moderate 
values of U. In the former case, the gap is induced by the 
external modulation and the shells consist of fractional and 
unit filled insulator domains surrounded by superfluid regions. 
In the later case, the gap is induced by the onsite interaction 
energy U and the shells consist of integer filled MI domains 
surrounded by superfluid regions. 




-40 



FIG. 1 1 : (color online) Density distribution for a system with A = 
1.1 and il/J — 0.005 as a function of the total number of atoms. The 
grid lines are at TVf'^ = 9, A^i = 32 and N^'^ = 42. 




40 



modes are populated, number fluctuations become only rele- 
vant at the edges. The insulator character of the atoms at these 
central sites and the reduced number fluctuations can be seen 
in Figim 

For N > Ni the trapping energy cost of placing an atom 
at the trap edge is higher than the energy needed to place the 
atom at the center For 1 < iV < N2^, atoms occupy the ex- 
tended modes of the second band and tend to spread over the 
different empty sites at the trap center. When the number of 
atoms N = iV^^, the first localized mode of the band is pop- 
ulated and the site at the trap center acquires filling factor one. 
As A^ is further increased, the width of the unit-filled central 
core grows. For N > A^^^ the density profile alternates from 
unit filled MI, superfluid, fractional MI and superfluid as one 
moves from the center towards the edge of the cloud. This 
shell structure in the HCB limit resembles in many aspects 



FIG. 12: (color online) Number fluctuations in the presence of trap 
for a system with A = 1.1 and Q/J — 0.005 as a function of the 
total number of atoms. The grid lines are at Ni^ — 9, Ni —32 and 

Ni'^ = 42. 

We now describe the first and second order interference 
patterns in the trapped system. The quasimomentum distri- 
bution is shown in Fig.^] Consistent with the behavior ob- 
served in the density and number fluctuations, for filling fac- 
tors N < N[^ and Ni < N < TV^-^ the quasimomentum 
distribution is sharply peaked at Q = and Qa = vr, signaling 
the superfluid character of the system. The insulating phases 
that appear for filling factor TVf ^ < A^ < iVi and N > A^;f ^ 
are signaled in the quasimomentum distribution by the drop of 
the peak intensities and the flattening of the quasimomentum 
distribution profile. 

Fig. ^Jshows the HBTI pattern for a different total number 
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FIG. 13: (color online) Quasimomentum distribution for a system 
with A/ J = 1.1 and Q./ J — 0.005 as a function of the filling factor. 
The grid lines are at iV/"^ = 9, iVi = 32 and Af|"^ = 42. 




FIG. 14: (color online) Noise interference pattern as a function of the 
filling factor or a system with X/J = 1.1 and Q/J = 0.005. The 
grid lines are at iVf'^ = 9, A''i = 32 and Ni''^ = 42. 



N, consistent with the superfluid properties of the system at 
these finings. Finally for N > iV^^^, a unit filled insulating 
domain at appears at the trap center This leads to a decrease 
in the amplitude of the central peak until it reaches a constant 
value. The peak at Qa = tt disappears as the unit filled core 
at the trap center grows. In the regime N > N.2^, Fig. ^| 
shows some additional small modulations of the central peak 
amplitude at certain fillings. The oscillations take place when 
a high energy mode of the first band, degenerated with another 
high energy mode of the second band, enters the Fermi sea. 

In analogy with the homogeneous case, the transition to in- 
sulating state can best be illustrated in the normalized inten- 
sity of the noise coiTelations, namely the visibility. In Fig. 1151 
we plot the visibility of the central and the superlattice peaks. 
As mentioned above, the normalized pattern is perhaps the 
best experimental observable as the normalization procedure 
filters some of the technical noise introduced in the measure- 
ment procedure. The normalization procedure maps the cen- 
tral peak maxima that appear before the formation of the insu- 
lator domains to minima. Furthermore, when normalized, the 
amplitude of the superlattice induced peak becomes compara- 
ble in intensity to that of the central peak. 



3.5 




of atoms. For N < N[^, only extended modes are occupied 
and the intensity of the central peak grows monotonically with 
N. Similarly the peak at Qa = it increases from the neg- 
ative value it takes for iV = 1 and becomes a peak with a 
positive intensity. As the localized modes enter the Fermi sea 
( A^ > Ni^) the rate of growth of the central peak slows down 
and this change in the slope generates a maximum in the inter- 
ference pattern. Note that in contrast to the strong reduction 
observed in the 51 = HCB case, in the presence of the trap 
only a decrease in the rate of growth is observed. The reason 
of this behavior is the fact that in the trapped case there is al- 
ways a superfluid component at the edges of the atomic cloud. 
For N[^ > N > Ni the interference peak at qa = tt is 
clearly visible, reflecting the staggered character of the phase. 
At A^ = A^i + 1, we see a rather sharp peak to dip transition 
at Qa — TT. Similar to the homogeneous system, the transi- 
tion signals the beginning of the population of an empty band. 
The dip becomes a peak when the number of atoms increases 
beyond a certain value. For Ni < N < N2^, one also ob- 
serves an increase in the growth rate of the central peak with 



FIG. 15: Visibility of the second-order interference peaks at Q = 
(top panel) and Qa — tt (bottom panel) as the total number of 
trapped atoms is varied. The grid lines are at Ni^ = 9, Ni — 32 

and Ni'^ = 42. 



VI. CONCLUSIONS 

In this paper, we have studied shot noise correlations for 
interacting bosons in the presence of an additional lattice. Al- 
though we have focused on a period-2 superlattice, this simple 
case exhibits all the key aspects of the phase diagram of more 
general superlattice potentials. 

One of the central results of this paper is that noise coiTe- 
lations may provide a practical tool to construct the phase di- 
agram of many body quantum systems exhibiting transitions 
from superfluid to fractional and integer Mott phases which 
also include the staggered insulating phase. Our study sug- 
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gests that the intensity of the second order Bragg peaks pro- 
vides a new order parameter to characterize these phase tran- 
sitions. Although our calculations in the soft core regime were 
done for systems with a reduced number of atoms and wells 
and further studies for larger systems are needed to confirm 
our conjecture, we believe that this result is important in the 
theory of quantum phase transitions. Furthermore, our study 
of the trapped system demonstrates the possible realization of 
these phases in the laboratory. 

Shot noise spectroscopy has already proven to be a useful 
experimental tool to identify Mott insulating states. We con- 
sider our analysis of the shot noise and its possible relevance 
in identifying various phase transitions in interacting bosons 
in superlattices is pertinent for atom optics experiments. Up 
to date all the experiments in super-lattice has been restricted 
to the Bose-Einstein-condensate regime. However, ongoing 
experimental efforts are trying to reach the strongly corre- 
lated regime in superlattices. The motivation of these efforts 
is not only to gain understanding of the many-body physics 
but also because atoms selectively loaded in super-lattices can 
have wider separation and stronger confinement, useful prop- 
erties for the implementations of lattice based quantum com- 
puting proposals. We hope that our study will stimulate fur- 
ther experimental probing of one dimensional bosonic sys- 
tems loaded in superlattices via HBTl. 



N . 



T'^f 



j-e/A 



H^f ^ _(^Hl^ +Hl^), 



H, 



e/ _ 



-2JV'2(aI +0.1 - Vi) 

-(2A + /i)ni + —hi{ni - 1), 



H^^ = -2JVi(4 + a2-V'2) 



(15) 
(16) 

(17) 



+ (2A - /i)n2 + —n2{n2 - !)• 

The above Hamiltonian can be diagonalized using a Fock state 
basis truncated below a certain occupation value. However, 
second order perturbation theory using ip as an expansion 
parameter can provide a good analytic approximation of the 
phase diagram. At zero order in ijj the Hamiltonian is diago- 
nal in the Fock basis. The integer occupation numbers rii, n2 
that minimized the energy are given by the conditions 



_, /i + 2A 
ni - 1 < — - — < ni, 

fi-2X 
n2- I < — — — < n2. 



We consider three different cases: 
• -2A < fi< min(t/ - 2A, 2A) 



(18) 
(19) 



VII. APPENDIX A 



In this case the ground state of the system is the Fock State 
JOIOIOI . . . ) and the unperturbed energy is Eg = — 2A — /i. 
The first order correction to the energy vanishes and to second 
order one gets: 



To study the phase diagram for the homogeneous system we 
use the well know decoupling approximation 12711 and substi- 
tute: 



alai 



= Viflj + i'jal ~ ^iV^j 



(13) 



Tij is the 



into the BH Hamiltonian. Here ipj — (dj) ? 
superfluid order parameter and rij is the expectation value of 
the number of particles on site i. We will assume the order 
parameter to be real. This substitution leads to : 



(2) _,ji J^'l Ji'l 



2J^l 



+ ViV'2 (20) 



^ V/^-ZA ^ + 2A ' M + 2A-[/ 



Minimizing the energy respect to ^1^2 yields the following 
equation: 



1 

2.J(ii+2\+U) 
{n+2X-U){fi+2\) 




(21) 



Eq. ( I21> has nontrivial solution only when the determinant of 
the matrix vanishes. Therefore the surface where the determi- 
nant vanishes determines the insulating phase: 



H'^f - -J^(V,+i+V,-i)(a]+%-V',) (14) 



U 



+ ^(2Acos[7rj] - ^)hj + —J2^:i(^J ~ ^ 



where we have also introduced the chemical potential /i. As- 
suming periodic boundary conditions, one can reduced the 
problem to a two mode system due the fact that the Hamil- 
tonian is exactly the same for all even and all odd sites. We 
define -02 = "021, "01 = V'2i+i, 02 = 0,21 and ai = 0,21+1, so 



4J2(^ + 2A + [/) = (/i + 2A - U){[f - 4A2), (22) 

where the bar denotes the dimensionless variables U = U/J, 
A = A/ J and /I = f^/J- In Fig- E]we show, with a crossed 
blue line, the solutions of Eq. (I22> in the jl vs. U plane, for 
A = 0.5, 4 and 8 (top right, bottom left panel and bottom right 
panels respectively). The region inside the blue loops corre- 
sponds to the 1/2 filled insulator. In the top left panel we also 
show with a blue dash-dotted line the critical value of t)c~ 
as a function of A. Uc~ corresponds to the smallest tJ in 
the lobe, below which the system is always a superfluid. Ex- 



actly at Uc 



v=\l2 



the upper and lower branches of the chemical 
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potential merge and the energy gap closes up. It is impor- 
tant to point out that in general mean field calculations do not 
accurately predict critical values but, in general, they capture 
very well the physics of the phase transition. 

• 2A < ^ < C/ - 2A 

In this case the ground state of the system is the unit filled state 
with exactly one atom per site. The zero order ground state 



energy is Eg — ~2fi. For this case we set -01 = -02 
second order the ground state energy is given by: 



■0. To 



4.) = 



SJ^V-^ 



sj^v-^ 



H~2X-U /i + 2A 

/4JV 4JV 
V^ + 2A /i-2A 



U 



4J0^ 



(23) 



Minimizing with respect to ip one gets an algebraic equation 
that determines the boundary between the superfluid and insu- 
lating phases. The solution is displayed in Fig|4]for A — 0.5, 4 
and 8. The critical value of C/^^^ is also shown in the top left 
panel with a dotted red line. Note the critical value increases 
with A. This increase is in agreement with the idea that for 
weak interactions disorder helps to delocalize the atoms. 

• U -2X< fi<2 min([/ - A, A) 

In this case the on-site repulsion is not large enough to avoid 
double occupancies and the ground state of the system is the 
state with two atoms in the low energy sites and zero in the 
others. This situation can not be described in the hard core 
regime and it is only present at moderate values of U. To zero 
order in J, the ground state energy is E, 
To second order it is given by 



4.) 



2,/,2 



4J^V^ 



2,/,2 



H-2\-U 

+4J-0102 



4J^ 
U + 2A 



4 


) _ 


-2/i 


- 4A + [/. 


4J' 


' / 2 

-A* 




8J2^2 


2A- 


M 


+ 2\-U 








(24) 



By first minimizing with respect to 0i_2, and finding the so- 
lutions for which the determinant vanishes one obtains the 
boundary for this phase. The solution is displayed in Fig|3 
for A — 0.5, 4 and 8 with a black line. At mean field level 
the minimum value of A required for the existence of this in- 
sulating phase is A = 3.75. For values of A > 3.75 the unit 
filled system is a superfluid for U < U^^^, a staggered insu- 
lator for JJ^^^ < U < U^2^ and a unit filled Mott insulator 



for tj > U^^^. In the top right panel, the lower and up- 
per branches of the solid black curve correspond to C/^j^^ and 
Uc-T^ respectively. 



VIII. APPENDIX B 

The super-lattice potential splits the main band into two 
different sub-bands. In this "multi"-band picture, it is sim- 
ple to understand the modification introduced by the trap. 
For excitation energies below the band width of the first sub- 
band: S(") - £;(") < El^lf^ EE 2\/j2^~A2 _ 2A, the modes 
are delocalized states that spread every other site symmet- 
rically around the potential minimum. In the semiclassical 



picture these modes only see the classical turning point at 

x^^il = ±ay^(£;(") + 2VA2 + j2)/0. For higher energies, 

^(») — £;(o) > E^^^f^, the eigenstates become localized on 
both sites of the potential. These modes see besides the clas- 
sical turning point, the Bragg turning point of the first band 

x^g'> = ±a\/(£;(") + 2A)/f}. 

At excitation energies higher than 4A (the energy gap be- 
tween the two subbands), the potential energy cost of local- 
izing a state at the edge is larger than the energy required 
to populate the high energy wells at the trap center. As a 
consequence, states with quantum number n > ^/4X/^i ap- 
pear centered again around the trap minima. These second 
group of extended states see the classical turning point at 

x^^l = ±ay'(£;(") - 4A + 2VW^~P)/n. Finally, for exci- 
tation energies larger than the band width of the second band 

eM - £;(«) > E^% 

— with 



^with^ the modes feel the Bragg 

.(2) 



turning point of the band x^g' ^ ±a^(£'(") - 2\)/n and 
become again localized at the edge. Sometimes a high energy 
mode in the second band becomes degenerated with some of 
the very high energy modes of the first band. This effect can 
be seen in Fig. [TO] 
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